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ABSTRACT 

Gravitational lensing of a background source by a foreground galaxy lens occasionally 
produces four images of the source. The cusp and the fold relations impose conditions 
on the ratios of magnifications of these four-image lenses. In this theoretical investiga- 
tion, we explore the sensitivity of these relations to the presence of substructure in the 
lens. Starting with a smooth lens potential, we add varying amounts of substructure, 
while keeping the source position fixed, and find that the fold relation is a more ro- 
bust indicator of substructure than the cusp relation for the images. This robustness 
is independent of the detailed spatial distribution of the substructure, as well as of 
the ellipticity of the lensing potential and the presence of external shear. 

Key words: theory — gravitational lensing — substructure 



1 INTRODUCTION 



Gravitational lensing describes the deflection of light rays 
from a background source caused due to the presence of a 
foreground concentration of mass (Refsdal 1964; Blandford 
& Narayan 1986; Schneider et al. 1992). The source can be 
a star, a galaxy, or a quasar, while the matter concentra- 
tion, the lens, can also be a star, an individual galaxy, or a 
cluster of galaxies. The geometry of space-time, defined by a 
particular cosmological model, is also an essential ingredient 
in gravitational lensing, as the strength of the deflection de- 
pends on the relative positions of the source, the lens, and 
the observer. If the lens is a very massive object, such as 
a galaxy or a cluster of galaxies, and if it is appropriately 
aligned with the background source, then it can produce 
multiple images of the source. This defines the 'strong lens- 
ing', regime. In this work, we study the strong lensing of a 
background source by a galaxy lens. The salient feature of 
gravitational lensing is that this mapping from the source 
plane to the image plane preserves surface brightness. How- 
ever, because the monochromatic flux of the source is usually 
unobservable, we can only measure the magnification ratios, 
or flux ratios, of the lensed images. Based on these flux ra- 
tios, we can then construct a model of the gravitational lens 
that will most accurately reproduce observed values. 

Unfortunately, this modeling has turned out to be very 
difficult in practice, despite stringent constraints. The pri- 
mary constraint on lens modeling comes from the magnifi- 
cation theorem, which states that for a given source posi- 
tion, the magnification of all images must sum to at least 
unity. One of the most important issues in gravitational lens- 
ing is the (growing) number of lenses whose flux ratios vio- 



late this magnification theorem. Some of the lenses with so- 
called 'anomalous flux ratios' are B1422+231, PG1115+080, 
B0712+472, B2045+265, and SDSS 0924+0219 (Patnaik et 
al. 1999, Chiba et al. 2005, Jackson et al. 1998, Fassnacht et 
al. 1999, Inada et al. 2003). These lenses, among others, ex- 
hibit flux ratios that cannot be fit with smooth lens models. 
It is reasonable to assume that these anomalous flux ratios 
may be due to limitations in our lens models. For example, 
our lens models may have inaccurate or too few parameters. 
However, as Mao & Schneider (1998) first demonstrated, for 
the case of B 1422+231, in which the discrepancy between 
observed and model-predicted flux ratios cannot be ascribed 
to an incorrect choice of parameterization, but rather to sub- 
structure in the lens. 

This suggestion is compelling, because it is in accord 
with predictions of the granularity of the mass distribution 
inferred from high resolution cosmological simulations in a 
Cold Dark Matter (CDM) dominated Universe. These sim- 
ulations have found copious amounts of substructure on all 
scales in the Universe, including galaxy scales (Diemand et 
al. 2005; Mathis et al. 2002; Moore et al. 1999). Metcalf 
& Madau (2001) quantified the effect of CDM substructure 
on flux ratios and found that 10 4 — 1O S M0 substructures 
near the Einstein radius can cause the flux ratios to deviate 
significantly from their model-predicted values. Soon after- 
wards, Dalai & Kochanek (2002) introduced a method by 
which to measure the abundance of substructure using lens- 
ing data, and concluded that substructure comprised ~ 2% 
of the mass interior to the Einstein radius of typical lens 
galaxies. 

Assuming, therefore, that substructures exist and that 
their effect is important, we would like to develop a robust 
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diagnostic that quantifies their presence. In the theoretical 
investigation presented in this paper, we keep the source 
position fixed and study the effect of substructure on the 
image configurations. Aided by previous work from Keeton 
et al. (2003, 2005), we have two diagnostics at our disposal, 
namely, examining the cusp and the fold relations (R C us P and 
-Rfdd hereafter), which we do in this paper. We compare the 
robustness of these two relations using the publicly available 
software Gravlens by Keeton (200f). We simulate a simple 
lens model with added substructure, with the aim of seeing 
how -Rcusp and i?f id gauge the presence of substructure as 
a function of its mass and position in the lens. 

This paper is organized as follows. In §2 we give a brief 
review of the magnification relations for folds and cusps, and 
describe the general properties of R CU sp and -Rfoid- In §3 we 
use -Rcusp and i?foid on a singular isothermal ellipsoid (SIE). 1 
These two sections also provide an overview of previous work 
by Keeton et al. In §4 we add substructure to our lens and 
investigate the values of R CU s P and 7?f id as the masses and 
positions of the substructure are varied. We also examine 
these relations when the ellipticity of the lens is varied and 
when external shear is added to the lens. We present our 
results and discuss its implications in §5. 



2 THE CUSP AND THE FOLD RELATIONS 

We begin with a brief review of the necessary lensing ter- 
minology. The lens equation for a gravitational lens sys- 
tem relates the impact parameter of the source's light ray 
on the lens plane L to the source's position on the source 
plane S, by taking into account the deflection of the ray by 
the lens mass. It can be written in dimensionless form as 
y = x — a(x), where y is the position of the source, x is the 
source position on the lens plane, and a(x) is the bending 
angle vector, which accounts for the deflection of the light 
ray. The lens equation can also be viewed 'in reverse', as a 
map fj : L — » S tracing the light ray backwards from the lens 
to the source plane. That is, we can view if as the assignment 
x i ► fj(x) = x — a(x). The inverse of the determinant of the 
Jacobian of this map, (det[Jac?|](3;)) _1 , for a lensed image 
at position x on the lens plane, gives the magnification of 
that image, and is conventionally referred to as the amplifi- 
cation matrix A(x). The magnification is a signed quantity: 
when it is negative, the image is called a saddle; when it is 
positive, it is called a minimum or a maximum. 2 From the 
definition of magnification that those positions x for which 
det[J&cfj\(x) = formally have an infinite magnification. 
The collection of all such points on the lens plane defines 
the critical curve. The corresponding collection of points on 

1 In this paper, as in Gravlens, the ellipticity is defined by e = 

1 — q, where q is the axis ratio (the ratio of the minor axis to the 
major axis). 

2 The terms 'saddle', 'minimum', and 'maximum' are standard 
in Morse theory. The number of minus signs appearing across 
the diagonal of the Hessian matrix of a Morse function gives the 
number of 'downhill' directions of that function. In particular, if 
there are no minus signs, then there are no downhill directions, 
and hence the function has a minimum at that point. See Petters 
et al. (2001) for a comprehensive treatment of the use of Morse 
theory in gravitational lensing. 



the source plane defines the caustic curve. Focusing on the 
source plane, the smooth portions of the caustic curve are 
called its folds, while the points where two abutting folds 
coincide are called its cusps. Typical examples of critical 
curves, caustic curves, folds, and cusps arc shown in Fig. 1. 
For a detailed treatment of these concepts, see Blandford & 
Narayan (1986), Schneider et al. (1992), and Petters et al. 
(2001). 

For folds and cusps, certain local relations between the 
magnifications of the multiple images are satisfied (Bland- 
ford & Narayan 1986; Schneider & Weiss 1992; Schneider 
et al. 1992; Petters et al. 2001; Gaudi & Petters 2002a,b; 
Keeton et al. 2003, 2005). For example, when the source lies 
asymptotically close to a cusp caustic (see Fig. 1), the three 
closely-spaced images (the so-called cusp triplet) should sat- 
isfy \ha \ — | Ms | + |mc| ~ 0, where /i; is the signed magnifica- 
tion of image i. For four-image lenses, Keeton et al. (2003) 
used this relation for the cusp triplet to define 

R _ \ha\ - \hb\ + Mcl F A -F B +Fc . . 

cusp " M + + Fa + Fb + Fc' [ > 

where Fi = F BTC \fii\ is the flux of image i if the source has 
flux Fare (we have essentially divided out by F SIC since it is 
unobservable, and are left with a dimensionless quantity), 
/is is the magnification of the middle image, and there is no 
need to specify whether it is a minimum or a saddle (for an 
SIE, it is a saddle if the source lies on the long axis of the 
caustic, and a minimum if the source lies on the short axis). 
The ideal cusp relation f? C us P — » is satisfied only when 
the source lies asymptotically close to the cusp caustic. To 
move beyond the asymptotic regime, Keeton et al. (2003) 
expanded the lens mapping in a Taylor series about the cusp 
to get 

Rcusp = -t- A CUS pd + ■ ■ • , (2) 

where d is the maximum separation between the three im- 
ages and the coefficient A CUS p is a function that depends on 
properties of the lens potential at the cusp point. In fact, 
Keeton et al. (2003) found that the properties of the lens 
that matter for A cusp are the ellipticity, higher-order mul- 
tipole modes, and external shear, whereas the radial mass 
profile is unimportant. Looking at eqn.(2), we see that as the 
source moves a small but finite distance from the cusp, i? cusp 
picks up a term to second order in d. Keeton et al. (2003) 
derived reliable upper bounds on R C us P , and concluded that 
these bounds would be violated only if the lens potential 
has significant structure on scales smaller than the distance 
between the images. 

A similar magnification relation holds when the source 
lies near a fold caustic (see Fig. 1). In this case, two of the 
four images lie closely-spaced together, straddling the criti- 
cal curve (the so-called fold image pair) . One of these images 
is a minimum and the other a saddle. When the source lies 
asymptotically close to the fold caustic, the fold image pair 
should satisfy |^i m in| — Msad ~ 0. For four-image lenses, 
Keeton et al. (2005) used this relation to define 

r> - l/^uiin| |/^sad| Fmin F sa d , , 

-Kfold = r j— j 1 — — — -= l-i) 

l^min T" |^sad| -^min T" -fsad 

Like its predecessor, the ideal fold relation i?f id — > is 
satisfied only when the source lies asymptotically close to 
the fold caustic. To move beyond the asymptotic regime, 
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Figure 1. The three types of image configurations: folds (top 
right), cusps (middle right), and crosses (bottom right), with the 
corresponding source position and caustic curve to the left. The 
lens potential used here is an SIE with ellipticity e = 0.5. 



these authors expanded the lens mapping in a Taylor series 
about the fold point to get 

RioU = + Afoiadi + ■ ■ ■ , (4) 

where di is now the distance between the two images in the 
fold image pair and the coefficient ^foid is a function that 
depends once again on those properties of the lens poten- 
tial listed above. As we approach a cusp in the asymptotic 
regime, |Af id| — ► oo, where the sign depends on whether the 
cusp is on the long or short axis of the caustic. For a cusp 
on the long axis, Atoid — * — oo because the middle image is 
a saddle, whereas for a cusp on the short axis, -Afoid — > +oo 
because the middle image is a minimum (remember that we 
are now looking at a cusp triplet, but using R{ a id instead of 
^cusp, and that -Rfoid is to be evaluated on minimum/saddle 
pairs of images) . This means that iZfoid breaks down asymp- 
totically close to a cusp, which is not too surprising because 
it is designed to be evaluated only on fold points. We men- 
tion this fact because in our calculations below we do indeed 
evaluate Rf \d for cusp points, but we are safe in doing so 
because our source sits a small but finite distance from the 
cusp, and is not asymptotically close. 

The more interesting feature of iifoid, as Keeton et al. 
(2005) discovered, is that the validity of the ideal fold rela- 
tion depends not just on how close the source is to the fold 
caustic, but also where precisely the source is along the caus- 
tic. This is reflected in the coefficient ^4foid, which takes on 
all values, both positive and negative, as one moves along the 
caustic (see Fig. 2). For this reason the authors introduced 
another variable, the distance 0I2 to the next nearest image 
(note that, as we are dealing with minimum/saddle image 
pairs only, neither di nor d^ will ever denote the distance 
between two saddles or two minima). They did so because 
both the source and the fold caustic are of course unob- 
servable, and therefore the value of Afoid, too, is unknown. 
Fortunately, the source's position is encoded in the image 
configuration: not in the separation d\ between the fold im- 
age pair, but rather in the distance d2 to the next nearest 
image. (To give an example: a source near a fold but not 



Figure 2. Caustic curve for an SIE with ellipticity e = 0.5. The 
numbers below are the values of Af id in the asymptotic fold 
relation -Rf D id = A{ i^di, evaluated at the points selected above. 
Note that j4f G id is monotonically decreasing from left to right, as 
one moves from the short axis to the long axis. As Appendix A of 
Keeton et al. (2005) discusses rigorously, Afoid is to be evaluated 
on the caustic, but describes the fold relation in the vicinity of 
the caustic. 



near a cusp will have di <C d-z ~ Re, whereas a source near 
a cusp will have d\ ~ di <C -Re, where Re is the Einstein 
radius.) 3 Hence the authors concluded that when deriving 
an upper bound on -Rfoid, beyond which one can infer the 
presence of small-scale structure, d-z must also be taken into 
account along with d\. The sensitivity of -Rfoid with respect 
to the source's position along the fold caustic will unfortu- 
nately make things difficult for us, and for this reason we 
will concentrate on cusp points only, not fold points; see §3 
below. 

Incidentally, in Figures 1 and 2, and in all of our cal- 
culations below, we use as our lens model an SIE with el- 
lipticity e = 0.5. Although we are restricting our theoretical 
analysis to an SIE, our results are nevertheless general, for 
the following reason. Keeton et al. (2005) generated a large 
ensemble of realistic lens potentials drawn from three differ- 
ent observational samples, with varying values of ellipticity, 
octupole modes, and external shear. For each lens poten- 
tial in their ensemble, they chose ~ 10 6 source positions, 
solved the lens equation using Gravlens, and then computed 
(di, di, Rfoid) for each minimum/saddle image pair. They 
then plotted R{ id on the (di,^) plane. This enabled them 
to extract the probability distribution of i?f D id for fixed val- 

3 The Einstein radius defines the radial scale for a lensing con- 
figuration. If a star lies exactly behind another one, then due to 
the symmetry, a ring-like image appears. This ring has an angu- 
lar radius Be and a linear radius Re = ^l$e> where is the 
angular diameter distance from the observer to the lens plane. 
For distant galaxies acting as lenses, the Einstein radius Re is of 
order 1 arcsec. 
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ues of di and di- Next, using an SIE with e = 0.5, they 
chose source positions such that di° ld = d'j usp = .46R c i n , 
and calculated Rfoid. They compared this value of Rfoid to 
that of each of the lens potentials in their ensemble, for the 
same values of di and d.2- Finally, for each minimum/saddle 
image pair, they found that the value of Rfoid calculated in 
the case of the SIE fell within the probability distribution of 
their ensemble, and thus concluded that each image pair of 
the SIE was consistent with lensing by a realistic potential. 
They repeated this procedure on the 23 known four-image 
lenses, and concluded that each case in which Rfoid lay out- 
side the probability distribution constituted strong evidence 
that small-scale structure was present for that particular 
lens. Given their result, we take as our 'archetypal smooth 
lens potential' an SIE with ellipticity e = 0.5. However, we 
also vary the ellipticity and external shear and evaluate their 
impact on our results. 



3 INVESTIGATING Rcusp AND R FO ld 

Prior to adding our point-mass substructures to the SIE, we 
calculate the values for R CUS p and Rfoid for a smooth lens 
model. In what follows, we assume that the source sits a 
very small but finite distance from the cusp or the fold. 

For a source near a fold caustic but not near a cusp 
(mathematically, this means that we pick a neighborhood 
around our source that does not contain a cusp point), the 
image configuration is given in the top panel of Fig. 1. Cal- 
culating i?foid for the fold image pair (A,B) that straddles 
the critical curve gives Rfoid ~ because, as stated in §2, 
A & B have roughly equal and opposite magnifications. Of 
course we can apply Rf Q id to any minimum/saddle image 
pair. Now, the two fold pair images A & B will have much 
higher magnifications than the other two images C & D. In 
fact, it can be as much as three orders of magnitude greater. 
So the fold relation for an image in the fold image pair and 
one of the other two images will converge to Rfoid ~ ±1- If 
the image in the fold image pair is a minimum, we get +1; 
if it is a saddle, we get — 1. 

The more interesting case is when a source lies near a 
cusp point, in which case we can use both relations R CU sp 
and -Rfoid; see the middle panel of Fig. 1. In this case the 
magnifications of the outer two images A & C will have the 
same sign and exactly the same magnitude, while the mag- 
nification of the middle image B will have the opposite sign 
and roughly twice the magnitude, with the sign depending 
on whether the source lies near a long axis cusp (B is a sad- 
dle) or near a short axis cusp (B is a minimum). Of course, 
as stated in §2, for the cusp triplet (A,B,C), R C us P ~ 0. We 
now apply -Rfoid to the cusp triplet. Pairing the middle image 
B with either of the outer images A or C gives .Rfoid w ±1/3. 
If B is a saddle, then A & C are minima, so we get —1/3; if 
B is a minimum, then A & C are saddles, and we get +1/3. 
The fourth image D, not part of the cusp triplet, will have a 
magnification much less than that of A,B, or C. It can be as 
much as three orders of magnitude less. So the fold relation 
gives -Rfoid ~ ±1 for the combination of any one of the cusp 
triplet images and D. The sign will depend on which cusp 
triplet image we use, and whether the source sits on the long 
or short axis: for a long axis cusp, (A,D) and (C,D) give +1, 
while (B,D) gives — 1. 



Of all these values, the important ones for our purposes 
are the relations R CU sp ~ and -Rfoid « ±1/3 for the cusp 
triplet: as we will show in the next section, substructure 
breaks the cusp triplet symmetry. The breaking of the sym- 
metry causes the two outer images to no longer have identi- 
cal magnifications (and thus the same value for -Rfoid). They 
are also no longer equidistant from the middle image. 



4 SIMULATING THE EFFECTS OF 
SUBSTRUCTURE 

Since our examination involves a simulated configurations 
rather than observational data, we can and choose to keep 
the source position fixed and investigate the effect of sub- 
structure on the image configurations. Specifically, we inves- 
tigate the change to both the image configurations and to 
the values of Rfoid and R CUS p when we distribute substruc- 
ture in the form of point-masses onto our archetypal smooth 
lens potential, an SIE with ellipticity e = 0.5. We consider 
both random and symmetric spatial distributions for 5-10 
point-masses. We examine cases when the spatial distribu- 
tion of these point-masses is (1) less than, (2) roughly equal 
to, and (3) greater than the distance between the images. 
We use Gravlens to solve the modified lens equation (SIE 
+ point-masses) for a source placed very close to either a 
fold or cusp. In each case, we gradually increase the mass of 
our point-masses from (the control case) to the 'cut-off' 
mass, which is the mass value beyond which we no longer 
have a four-image lens. Finally, we calculate Rfoid and R CUS p 
for these granular lenses. 

When considering sources near a cusp, we concentrate 
on long axis cusps because they satisfy the ideal fold rela- 
tion more readily than short axis cusps (Keeton et al. 2005). 
We also ignore cross configurations (see the bottom panel of 
Fig. 1) because in such a case di ~ d-z ~ Re, and scales on 
the order of the Einstein radius are no longer 'small-scale'. 
(Recall from §2 above that a violation of the cusp or fold 
relation implies the presence of significant substructure on 
scales smaller than the distance between the images. There- 
fore calculating Rfoid for widely separated images does not 
really tell us about small-scale structure, because we are no 
longer restricted to short length scales.) 

Of course, the shortest length scale (i.e. the smallest 
possible value for di) is the distance between the fold image 
pair in a fold image configuration (the two closely-spaced 
images in the top panel of Fig. 1). However, we encounter a 
problem when investigating sources near a fold but not near 
a cusp. In general, changing the lens potential by adding 
substructure changes the shape of the caustic curve, and 
the nature of the change depends on both the spatial dis- 
tribution of the substructure, and more importantly, their 
masses. To give an example: setting the mass of our SIE to 1 
(this sets Re = 1 in the code units) , the spatial distribution 
of five point-masses in the manner shown in Fig. 3, with 
masses ranging from < rrii < 0.0190, deforms the caustic 
curve considerably, but still produces four-image lenses. For 
masses 0.0190 <m, < 0.0259, the lens can also produce five 
images, though the fifth image is quite faint: its magnifica- 
tion is \fi\ ~ 10~ 5 , while for any one of the cusp triplet 
images is ~ 10 2 . Finally, for masses mi > 0.0259, the lens 
can produce six images. Strictly speaking, the cusp and the 
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fold relations are defined only for four-image lenses, so we 
will not examine -R CU s P and R{ id in the five- and six-image 
regimes. 

We would like to come as close as we can to satisfy- 
ing the ideal fold and cusp relations, we want to place our 
source as close to the caustic curve as Gravlens allows. Now, 
aside from deforming certain regions of the caustic curve, the 
addition of substructure also tends to narrow the diamond- 
shaped region of the curve, so that a source placed very close 
to a fold or cusp may actually end up falling outside the 
diamond-shaped region of the new caustic curve. A source 
outside this region, of course, will not produce fold or cusp 
configurations, so this means that we are forced to move 
our source back inside the diamond-shaped region of our 
new caustic curve in order to produce four-image configura- 
tions. We do encounter the problem of trying to place our 
new source in the 'same place' with respect to the new fold 
caustic as our old source was with respect to the old fold 
caustic. Given the sensitivity of iifoid to where our source 
lies along the fold (see §2 and Fig. 2) , this is difficult in prac- 
tice. When we move our source back inside the caustic curve, 
the value of i?f id will change, but we cannot say with con- 
fidence whether this change in 7?foid is due to the addition 
of substructure, or simply because we displaced our source 
by a small amount. For this reason, we ignore sources lying 
near a fold. 



Fortunately, we do not encounter this problem when 
we place our source near a cusp, because the cusp is a much 
smaller region than the fold caustic. We cannot 'move along' 
the cusp in the same way that we can move along the fold 
because the cusp is essentially a tiny wedge. Therefore we 
can safely move our source back inside the caustic curve and 
place it back into the wedge of the cusp. We still have to be 
careful, however, because the addition of substructure tends 
to displace the long axis cusp slightly off the axis (so that 
the long axis is no longer aligned with the coordinate axis), 
making it more difficult to locate the cusp accurately. Doing 
so gives us our most interesting result: when we add sub- 
structure onto our SIE, and consider a source near our new 
long axis cusp, we no longer have the symmetry in our cusp 
triplet. The cusp triplet tends to be displaced to one side, so 
that the two outer images are no longer identical. In fact, one 
of them is closer to the middle image. This is illustrated in 
the three panels shown in Fig. 3. At this point an important 
issue arises: because the cusp point is more difficult to locate 
when substructure displaces the long axis, one may wonder 
whether the skewed cusp triplet shown in Fig. 3 is caused 
merely by the possibility that we have missed the cusp point 
and instead placed our source against a fold caustic. If so, 
then we would be able to reproduce the same skewed cusp 
triplet in the case of an SIE without substructure, simply by 
displacing our source off the cusp by a small amount. But 
in fact this is not the case: a source displaced slightly off 
the cusp for an SIE without substructure produces a very 
tight fold configuration (as in the top panel of Fig. 1), one 
that looks altogether different from the skewed cusp triplet 
shown in Fig. 3. Thus it is the substructure that breaks the 
cusp triplet symmetry. 




Figure 3. The breaking of the cusp triplet symmetry for a distri- 
bution of five point-masses onto an SIE with cllipticity e = 0.5. In 
each panel, the diagram on the left shows the caustic curve and 
source position while that on the right shows the corresponding 
critical curve and cusp configuration. For each caustic curve we 
have moved our source as close to the cusp as Gravlens allows. In 
each case, note the deformations of both the caustic and critical 
curves due to the point-masses, and compare these with the cusp 
configuration in Fig. 1. The triangles are the point-masses and 
the squares are the images. In each panel, the point-masses have 
the same mass. In the first panel, that mass is m = 0.0185 (in 
units where = 1); in the second panel, m = 0.0255; in the 
third, m = 0.0258. Note the fifth image near the origin in the 
second and third rows. For m > 0.0259, we get six images. 

5 RESULTS AND CONCLUSIONS 

We now examine the sensitivity of the cusp and the fold 
relations to the spatial distribution of substructure shown 
in Fig. 3. We set the mass of our SIE to I (Re = 1) and 
vary the substructure masses from (the control case) to 
Hij = 0.019. For a typical galaxy lens this corresponds to 
adding substructure in the mass range of 10 5 - 1O 7 M . Note 
that for m,i > 0.019 we no longer have a four-image lens. 

5.1 Dependence on Substructure Mass 

First of all, for the control case of the smooth poten- 
tial and no substructure, we find R CU s P = 0.002 and 
iifoid = (-0.331, -0.331), where J? fo i d = -0.331 for both 
the left/middle image pair and the right/middle image pair. 



6 Aazami & Natarajan 



As discussed in §3, we expect a value close to for -R CU s P 
because we have placed our source as close to the cusp as 
Gravlens allows (see Fig. 1). However, our source still sits 
a finite distance from the cusp, not asymptotically close, 
so we should not expect to satisfy the ideal cusp relation 
Rcusp — > exactly; hence R C u Sp = 0.002 is acceptably close 
to zero. As for Rfoid = —0.331, this value, too, is what we ex- 
pect: both the leftmost and rightmost images have the same 
magnification, which is half the magnitude as that of the 
middle image with the opposite sign. In this case our source 
sits near a long axis cusp, so the middle image is a saddle 
and -Rfoid is negative. As discussed in §3, we expect a value 
of -Rfoid ~ —0.333, so a value of —0.331 is close enough. 

Now we examine the breaking of the cusp triplet sym- 
metry shown in Fig. 3, for three mass values in the four- 
image regime: 

= .002 and R iold = (-.331, -.331) 
= .005 and R fo i d = (-.326, -.332) 
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m = 


.000 
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m = 


.001 


^cusp 
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m = 


.005 


^cusp 


• 


m — 


.010 


-^cusp 



flcusp = .006 and R io id = (-.317, -.340) 
.006 and R fold = (-.305, -.352) 

The notation R fo i d = (-0.317, -0.340) implies that the 
left/middle image pair gives Rioid = —0.317, while the 
right/middle image pair gives Rioid = —0.340; likewise for 
the others. The change in Rfoid manifestly reveals the break- 
ing of the cusp triplet symmetry; the change in R C us P , how- 
ever, does not, as shown in the top panel of Fig. 4. In Fig. 4 
we show the effect of varying the mass of the substructure 
on i?cusp and .Rfoid- It is clear that, within the four-image 
regime, Rfoid readily reflects the change in the mass of the 
substructure, while R C us P stays roughly constant. 

Keeton et al. (2005) demonstrated that, for observed 
cusp lenses, the fold relation indicated the same flux ra- 
tio anomalies as the cusp relation, with one exception: in 
B2045+265, the possible presence of octupole modes pre- 
cluded the authors from declaring that the fold relation was 
violated, whereas the cusp relation was very clearly violated. 
However, because the fold relation is defined for pairs, it 
gave the additional information of indicating which partic- 
ular image in the cusp triplet was the one most affected by 
substructure, a distinction that could not be made with the 
cusp relation. What we have shown in our analytic argu- 
ment here is that, for an SIE with e — 0.5, the fold relation 
is a more reliable indicator of substructure than the cusp 
relation. 



5.2 Sensitivity to External Shear 

Both i?cusp and -Rfoid are expected to be sensitive to the pres- 
ence of external shear. We therefore add external shear with 
an amplitude of 7 ~ 5% to the configurations. With this 
modification, we examine the breaking of the cusp triplet 
symmetry for the same substructure mass values listed in 
§4.1, and find: 



m = .000 : -Rcusp 
m = .001 : i?cusp 
m = .005 : R C u Sp 
m = .010 : 



.005 and Rfoid 
.006 and Rfoid 
.006 and Rfoid 



(- 



Rcusp = -006 and R fo i d 



329, -.329) 
.326, -.330) 
= (-.318, -.338) 
= (-.311, -.346) 



(- 



Once again, -Rfoid manifestly reveals the breaking of the cusp 
triplet symmetry, whereas R C u Sp does not. In fact, as Fig. 4 
shows, -Rcusp eventually falls below the control case value of 
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Figure 4. i? C usp and Rf Q \d as a function of the mass m of the 
substructure for the e = 0.5 case, for the spatial distribution given 
in Fig. 3. Left column: Without shear. For < m < 0.0185 (in 
units where = 1), we get four images for a source as near the 
long axis cusp as Gravlens allows. Right column: External shear 
with an amplitude of 7 ~ 5%. In this case we have a four-image 
regime for < m < 0.020. The control case (no substructure) 
is shown by the dashed line. In the bottom panel, the black and 
white data points denote the left / middle and right / middle image 
pairs in the cusp triplet, respectively. i?foi d reflects the change in 
the mass of the substructure better than -Rcusp- 



Rcusp = 0.005. On the other hand, -Rfoid remains correlated 
with the mass of the substructure throughout the four-image 
regime, even with external shear. Finally, note that the four- 
image regime in this case is produced when m; < 0.02. 



5.3 Sensitivity to Ellipticity 

The ellipticity of our SIE affects the values of R C u Sp and 
Rfoid ■ We find that for ellipticities lower than our archetypal 
value of e = 0.5, Rcusp remains unresponsive to the substruc- 
ture, whereas Rfoi d is still sensitive to the mass of substruc- 
ture, though not as robustly as for the case with e = 0.5. The 
left-hand panels of Fig. 5 demonstrates this explicitly for an 
ellipticity e = 0.25. For ellipticities higher than our archety- 
pal value, Rcusp displays erratic fluctuations, whereas -Rfoid 
remains well correlated with the mass of the substructure, 
as shown in the right-hand panels of Fig. 5. 



5.4 Concluding remarks 

In this theoretical investigation, we attempt to develop a ro- 
bust diagnostic that quantifies the presence of substructure 
in lens galaxies. To this end, we simulated a simple lens po- 
tential with added substructure, with the aim of seeing how 
the cusp and the fold relations, calculated while keeping the 
source position fixed, respond to the presence of the sub- 
structure as we vary its mass and position within the lens. 
We took as our lens model an SIE and used point-masses as 
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Figure 5. Left column: i? CUS p and Rfoid as in Fig. 4, but now 
for a lens with ellipticity e = 0.25. We have a four-image regime 
for < m < 0.009. We ignore the region < m < 0.002 be- 
cause it is difficult to cleanly place the source near the cusp 
point. Right column: The relations for a lens ellipticity e = 0.75. 
In this case we have a three-image regime (a cusp triplet, no 
fourth image) for < m < 0.035, and four-image regime for 
0.035 < m < 0.100. Note that although -R CU sp fluctuates errat- 
ically, in this case the right /middle image pair of R{ a \d is not 
reliable, either. The most reliable indicator of the substructure 
masses here is the left / middle image pair of if,f id . 
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our substructure with masses ranging from < m; < 0.019 
(Re = 1). When we varied the ellipticity of our SIE, we 
found that for low ellipticities the cusp relation was unre- 
sponsive to the substructure, whereas for high ellipticities 
it became erratic. The fold relation, on the other hand, re- 
mained well correlated with the mass of the substructure. 
Considering the effect of external shear gave the same result: 
the fold relation remained well correlated with the substruc- 
ture, whereas the cusp relation did not. 

We considered random distributions of point-masses as 
substructure and cases when their spatial distribution was 
less than, roughly equal to, and greater than the distance 
between the images. Overall, we found that the addition 
of substructure breaks the symmetry of the cusp triplet, 
and that the fold relation responds more accurately to the 
change in mass of the substructure than the cusp relation. 
We conclude, therefore, that the fold relation is the more 
robust diagnostic of substructure. 

In order to apply this technique to real data, we will 
need to use the observed image positions and the distances 
between them, and then use Monte Carlo methods to con- 
strain the corresponding source positions. While this is be- 
yond the scope of this work, we pursue it in a follow-up 
paper. 



